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Abstract 

Polynomial sequences Pn{x) of binomial type are a principal tool in 
the umbral calculus of enumerative combinatorics. We express Pn{x) 
as a path integral in the “phase space” N x [—7r,7r]. The Hamilto¬ 
nian is h{(f>) = it produces a Schrodinger type 

equation for Pn{x). This establishes a bridge between enumerative 
combinatorics and quantum field theory. It also provides an algorithm 
for parallel quantum computations. 
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Under the inspiring guidance of Feynman, a short¬ 
hand way of expressing—and of thinking about—these 
quantities have been developed. 

Lewis H. Ryder Chap. 5]. 


1 Introduction on Polynomial Sequences of 
Binomial Type 

The umbral calculus |^, ^ in enumerative combinatorics uses poly¬ 

nomial sequences of binomial type as a principal ingredient. A polynomial 
sequence Pn(x) is said to be of binomial type []^ p. 102] if Pq{x) = 1 and if 
it satisfies the binomial identity, 


k=0 


Pnix + y) = ^ Pk{x)Pn-k{y) 
for all n G N and x, t/ G M. 


( 1 ) 


Example 1.1 Among many examples [^, § 5] of polynomial sequences of 
binomial type we mention only very few: 

1. The power monomials Pn{x) = x'^. 

2. The rising factorial sequence Pn{x) = a;(a; -f 1) ■ ■ ■ (a; -b n — 1). 

3. The falling factorial sequence Pn{x) = x{x — 1) ■ ■ ■ (a; — u -b 1). 

4. The Abel polynomials An{x) = x{x — an)"'~^. 

5. The Laguerre polynomials 



Using delta functionals one can show that 
nomial sequence of binomial type Pn(x): 

Pn{0) = <^"0 = I q’ ^ ^ q’ (3) 

For a polynomial sequence Pn{x) we could consider its coefficients Ony in hs 
decomposition over power monomials x^: 

n 

Pn{x) = y^^OnyX^- 
k=0 


18, p. 106] that for any poly- 
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The coefficients Uny are an example of connecting constants |1l8| , § 6] between 
sequences Pn{x) and x”. The dehnition (|) of polynomial sequences of bi¬ 
nomial type can be expressed in the term of coefficients an,k as follows 
Prop. 4.3]: 

t + j\ _ I n 

i j ~ 2^ [k 


18, 




(4) 


Therefore these coefficients are highly interdependent. Particularly from (^) 
it follows that 0^,0 = ^n,o- Moreover one can observe from that the 
knowledge of the sequence c„ = a„_i, n = 0,1,... allows to reconstruct all 
an,k by a recursion. Sequence c„ = is known as cumulants associated 
with polynomials Pn{x) and the related probability measure, see e.g. |]T^ for 
a survey. One may state the following 

Problem 1.2 Give a direct expression for a polynomial sequence Pn{x) of 
binomial type from the sequence of its cumulants c„ = an,i = P^(0), n = 

0 , 1 ,.... 

There are several well known solutions to that problem. The umbral 
calculus suggests the following procedure[| (see 
sition). Let us dehne the formal power series: 


for an interesting expo- 


/W = 

71=1 


nl 


(5) 


then the function is a generating function for polynomials Pn{x): 


00 






( 6 ) 


lem 1.2 


The main goal of this paper is to present a different solution to the Prob- 
It is probably not the best from the computation point of view unless 


computations are made in aprallel by a quantum computer, see Remark ^ 


We will express Pn(x) as a path integral (|23|) in the “phase space” N x [—tt, tt]. 
The Fourier transform h(0) = of the sequence Cn/n\ = p'^{0)/n\ 

plays the role of a Hamiltonian. Thus we establish a connection between 
enumerative combinatorics and quantum field theory. This connection could 
be useful in both areas. It is interesting to note that quantum mechani¬ 
cal commutation relations PQ — QP = I was applied to umbral calculus 
by J. CiGLER H. Existence of a “yet-to-be-discovered stochastic process” 


am grateful to the first anonymous referee for its indication. 
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connected with Lagnerre polynomials was conjectnred by S.M. Roman and 
G.-C. Rota |]T^. Note also that the approach to the nmbral calcnlus de- 

^ can be a prototype for a general constrnction of 


veloped in 

coherent states {wavelets) in Banach spaces 0 


2 Preliminaries on Path Integrals 


Let a physical system has the configuration space Q. This means that we can 
label resnlts of onr measurements performed on the system at any time to eM. 
by points of Q. For a quantum system dynamics the principal quantity is the 
propagator [|, § 2.2], |2^, § 5.1] K{q 2 ,t 2 ',(li,ti )—a complex valued function 
dehned on Q x M x Q x M. It is a probability amplitude for a transition 
Qi <l 2 from a state qi at time ti to q2 at time t2- The density of probability 
of this transition is postulated to be 


P{q2,t2;qi,ti) = \K{q2,t2;qi,ti)[ 


The fundamental assumption about the quantum world is the absence 
of trajectories for a system’s evolution through the conhgurational space 
Q: the system at any time ti could be found at any point q^. Then the 
transition amplitude qi —g 2 is a result of all possible transitions amplitudes 
qi ^ qi ^ <?2 integrated over Q. 

Mathematically this can expressed in the term of propagator K as follows. 
Let us £x any U, ti < ti < t 2 . There exists a measure dq on Q such that 


K{q2,t2-,qi,ti) 


K{q 2 , ty qi, ti)K{qi, tp, qi,ti) dqi. 


(7) 


This identity could be thought as a recurrence formula for the propagator 

K. 

Another assumption on propagator is its skew symmetry 

P{q2,t2]qi,ti) = K{qifiyq2,t2), ( 8 ) 

it ensures that all transitions qi —> q 2 are done by unitary operators. 

It is useful to consider K{q 2 , ^ 2 ; <?i, R) as a kernel for the associated integral 
operator K : ^ 2 ( 0 , dq) —> ^ 3 ( 0 , dq) defined as 

[770](g,t)= / K{q,t;qi,ti)(j){qi,ti)dqi, (9) 


with ti is hxed. Operator K is defined on functions (j){q,t) G L2{Q,dq) such 
that the right-hand side of (j^) does not depend on ti. Then (|^ is read as 
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= K and implies K* = K. By other word K is an orthoprojection[| 
from L 2 {Q,dq) to its proper snbspace Imi^. 

Another conseqnence of (|^): fnnctions (j){g 2 ,t 2 ) which belongs to image 
of operator K are fixed points for K-. 

(f){q^t) = [K(t)\{q,t) = ! K{q,t;qi,ti)(t){qi,ti)dqi. (10) 

Jq 


These are wave functions and describe actnal states of the system. It follows 
from that all fnnctions n{q,t) = K{q,t] q'o,to) for arbitrary hxed qo G Q, 
to G M are wave fnnctions too. The fnnction \(j){q,t)f gives the probability 
distribntion over Q to observe the system in a point q at time t if it was in 
qo at time to. 

R. Feynman developing ideas of A. Einstein, Smoluhovski and 
P.A.M. Dirac proposed an expression for the propagator via the “integral 
over all possible paths”. In a different form it can be written as follows |^ , 
(5.13)1 




VqVp 
~~h 


■ exp 



dt {pq - H{p,q)) 


( 11 ) 


Here H{p, q) is the classical Hamiltonian [|I], Chap. 9] of the system—a fnnc¬ 
tion defined over the phase space. In the simplest case the phase space is the 
Cartesian prodnct of the configuration space Q and the space of momenta 
dual to it. Points in phase space have a double set of coordinates (g, p), where 
q is points of the configuration space Q and p is a coordinate in the space 
of momenta. The inner integral in (|^ is taken over a path in the phase 
space. The outer integral is taken over “all possible paths with respect to a 
measure VqVp on paths in the phase space”. We construct an example of 

similar to (ED in the next section. 

one can obtain an impression 


22 


path integral 

From reading of physics textbooks 
that the measure VqVp could be well dehned from an unspecified compli¬ 
cated mathematical construction^. But this is not true up to now: there is 
no a general mathematically rigorous definition of the measure in the path 

^For a model of relativistic quantisation generated by two orthoprojections see 
^The editor preface to the Russian translation of [D tells: “There is no a rigorous 
definition of the path integral in the book of Feynman and Hibbs... However it is not 
very important for a physicist in most cases; he only need an assurance that a rigorous 
proof could be obtained”. The recent textbook which were already reprinted 9 times 
in two editions refereed for a rigorous path integral construction to the old paper |^, 
however it is known that this approach is not complete, see e.g. the review of this article 
MR H 17,I26Ic as well as Q and ||^, Notes X.II]. 
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space. However there are different approaches of regularisation of path inte¬ 
grals, see | 2 ^ for recent accounts. 

It is known that path integrals give fundamental solutions of evolution¬ 
ary equations [^. Particularly wave functions (|I0[) are solutions to the 
Schrodinger equation |^, § 5.2], i § 4.1] 


ih 


d(piq,t) 

dt 




( 12 ) 


where H{q,p) is an operator which is the quantum Hamiltonian. It was 
mentioned that functions K{q, t) = K{q, t] qo, to) are wave functions and thus 
are solutions to ( 0 )- In fact it follows from (|^) that they are fundamental 
solutions, i.e. any other solution is a linear combination {superposition in 
physical language) of those. Therefore the propagator K{q 2 , ty, qt, U) contains 
the complete information on dynamics of a quantum system defined by the 
equation ([T2|). 


3 Polynomial Sequences from Path Integrals 

Let the conhguration space Q be a semigroup with an operation denoted by 
-|- and let a propagator K{q 2 ,t 2 ,qi,ti) be homogeneous in time and space, 

i.e. 

^{q2 + 9) ^2 + qi + <?; H + ^) = K{q2i ^ 2 ; qi, ti), (13) 

for all q,qi,q 2 ^ Q, t,ti,t 2 G M. Such a propagator corresponds to a free 
time-shift invariant system. Then K{q 2 ,t 2 ,qi,ti) is completely defined by 
values of the function S{q,t) such that S{q2 — qi,t2 — ti) = K{q2,t2,qi,ti). 
From (|^) we can deduce a “recurrence” identity for the function S{q,t) [^: 

S{q,ti + t2) = [ S{qi,ti)S{q - qi,t2) dqi. (14) 

Jq 


An equation of the form (H) was taken in [Hj as the dehnition of tokens 
between two cancellative semigroups Q and M’*'. It allows to express a time 
shift by the integral of all space shifts over Q. Tokens are closest relatives to 
sectional coefficients considered by Henle [0. A polynomial sequence Pn{x) 
of binomial type is another example of token O: the sequence qn{x) = 
Pn{x)/n\ is a token between cancellative semigroup N and M, i.e. from (|l|) 
follows a realization of (|^ in the form 


qn{x + y) = ^ qk{y)qn-k{x), 


k=0 


(15) 
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i.e. in this case the conhguration space Q is the discrete set of nonnegative 
numbers and the integration reduces to summation. 

If polynomial sequence qn{x) represents a particular propagator for a 
conhguration space Q = N then one may wish to obtain a path integral 
representation in the form (|Tl|) for it. We do it in a way parallel to the path 
integral deduction in |2^, § 5.1]. 

Let X G M be hxed and iV G N be an arbitrary positive integer. Then 
repetitive use of ([T5|) allows to write 

^1+^2 H- \-kN=n 

The hrst two terms in the Taylor expansion of j = 1,..., iV are: 

(iv) ^ 

where the Kronecker delta as the hrst term come from (j^). We use the 
Pontrjagin duality between Z and the unit disk T = [—vr, tt] to construct an 
integral resolution of the Kronecker delta: 

TT 

4o = ^/ (18) 

— 77 


It is also convenient to introduce the Fourier transform h{p) for the numerical 
sequence g(.(0) from the same Pontrjagin duality: 


h(p) = 

k=0 


q'M 


1 


j e-^P^h{p) dp. 


— 77 


(19) 


Then we can transform the Taylor expansion (^) as follows: 

®(v) “ + °(v) 

77 77 

= tLf f e-‘r'‘h(p)dp^ + o(^') 


271 


271 


1 


e-‘-‘(l + Mp)^)<ip+o(^) 

— 77 
77 

— / exp (-ipi + dp + o(^) , 


( 20 ) 
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where the last transformation is based on the estimation 

<!xp(Mp)^) = (i + Mp)^)+£>(^) 


Substituting N copies of (poD to (|16|) we obtain 


Inix) = ( 1 ^) (- 1 ) ' " * 


ki-\-k2-\ - ^ 

1 


nJ 


X \ 

n) 


E-'E 

fciH- \-kff=n 


(2i) 


N 


exp (^-ipiki + dpi 


X ■ ■ ■ X / exp 


-ipNkN + Hpn)^] dpN + o( ^ ) 


N 


\nJ 


7T TV 


E-'E 

ki-\ - \-kj\[=n 


( 2 vr: 


iN 


exp (^-ipiki + 


— TV —TV 


E-'E 

fciH- \-kj^=n 


X ■ ■ ■ X exp (-ipNkN + h{pN)^^ dpi... dp^ + 

exp -ipih + \ ‘I'Pi -‘I'Pn 

( 21 ) 


(2i) 


TV TV 


— TV —TV 


.1 = 1 


+0 


Now we can observe that the expression XlILi {—Ipiki + h{pi)^^ looks like an 
integral sum for the function —ip{t)k'{t) + h{p{t)) on the interval t G [0,a:]. 

At this point we say together with physicists § 5.1] magic words and 
(^) “could be symbolically written for N —>■ oo” as 

X 

qn{x) = / VkVp exp / {-ip{t)k\t) + h{p{t))) dt, (22) 


where the hrst integration is taken over all possible paths k{t) : [0,a;] —>■ 
such that fc(0) = 0 and k{x) = and the path p{t) : [0,a:] —[—vr, tt] is 
unrestricted. It is enough to consider only paths k{t) with monotonic grow— 
other paths make the zero contributions. Indeed the identity qi{x) = 0 for 
I < 0 (made by an agreement) implies that contribution of all paths with 
k'{t) < 0 at some point t vanish. Here (and in the inner integral of (p^)) 
k'{t) means the derivative of the path k{t) in the distributional sense, i.e. it 
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is the Dirac delta function times ^ in the points where k{t) jumps from one 
integer (mod value k(t — 0) = ^ to another k(t + 0) = 

In fact there is no any magic in the above path integral. Due to the 
property qn{x) = 0 for n < 0 paths k{t) take values on the finite set fco, 
fci,..., kn, otherwise they make the zero contribution to the path integral as 
was mentioned above. Therefore the set of paths k{t) is countable. Under 
such a condition the situation is equivalent to a finite Markov process with 
a well-dehned measure on the path space according to Kolmogorov theorem, 
see for example [^, § 1] or P5| , Th. 2.1]. Moreover at the beginning of the 
next section we show how to calculate the path integral (p2| ) as an ordinary 
integral. 


Now we can give a solution to Problem |1.2| using the relation p„ 
n\qn{x) between Pn{x) and qn{x): 


[X = 


Pn{x) =n\ VkVp exp / {—ipk' + h{p)) dt, h{p) = 


fc =0 


ypk 

Id' 


(23) 


Note that a comparison of our Hamiltonian h{p) with the formal power series 
f{t) (^ made out of cumulants shows that h{p) = /(e*^). 


Example 3.1 Here is a list of Hamiltonians producing the polynomial se¬ 
quences from Example 1.1. 


1. The power monomials Pn{x) = x'^ are produced by h{p) = e*^. 

2. The rising factorial sequence Pn(x) = x{x + l) ■ ■ ■ {x + n — 1) is produced 
by the Heaviside function h{p) = x{p) = 

3. The falling factorial sequence Pn{x) = x{x — l) ■ ■ ■ {x—n+1) is produced 
by the shifted Heaviside function 


Kp) = x(p - tt) = ^ 

k=l 


k 


CO 


E 


^i{p-TT)k 

k 


4. The Abel polynomials An{x) = x{x — an)'^ ^ are produced by h{p) = 

^cxD ^ipk 

2^k=l k\ ^ ■ 

5. The Laguerre polynomials (^) are generated by the delta function 

Kp) = Kp) = 
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4 Some Applications 

Let us start from the question^ of calculational applicability of the path inte¬ 
gral formula (^31). We transform this path integral into an ordinary integra¬ 
tion. Because we could not make any seriously restricting assumptions about 
the function h{p) the following arguments have an algebraical-combinatorial 
nature rather than analytic. They are similar to heuristic manipulations with 
formal power series which are used in a deduction of formula 

For a convenience we scale by ^ the parameter t in the inner integral 

X 1 

j {-ip{t)k'{t) + h{p{t))) dt = j i-ipiiti)k[{ti) + xh{pi{ti))) dti, 

0 0 

where on the right-hand side new paths are fci(ti) : [0,1] —>■ N with the 
endpoint /ci(0) = 0 and ki{x) = n, and the path pi(ti) : [0,1 ] —>• [—tt, vr] is 
unrestricted. We continue renaming variables of integration ti, ki, and pi by 
f, k and p respectively. 

As was mentioned before k{t) is an increasing step functions uniquely 
dehned by a collection of n points tj G [0,1] (some of them could coincide) 
where it jumps by 1. Thus k'{t) = — tj) and consequently: 

\ n 

/ —ip{t)k'it) dt = — 'E p{tj). 


Thereafter we could transform (E3f) as follows: 


Pnix) = n\ VkVpexp {—ip{t)k'{t)+xh{p{t)))dt 


= n\ I VkVp exp j —+ j xh{p{t))dt ) (24) 


Now we use the rotational symmetry of the unit circle [— tt, tt] and that 
paths p{t) are completely unrestricted, therefore the whole set P of such paths 
is also rotational invariant. Thus the integration over that set P acts just like 
an averaging over the unit circle [— vr, tt] . On the other hand a random path 
k{t) is dehned by a random collection of n points tj G [0,1] consequently the 

H am grateful to the second anonymous referee who suggested this question. 
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integration of p{tj) for any j over all paths k{t) will produce n independent 
equal uniformly distributed random variables on the unit circle. Thus we 
could hnally express our path integral (^) as a simple integration: 


= n\ VkVp exp [-zJ2pit,)+ xh{p{t))dt 


i=i 


= nl 


J exp i dp 


TV 

n\ j dp. 

— TV 


(25) 


If we recall the connection h{p) = /(e®^) between our Hamiltonian h{p) 
and the generating function of cumulants f(t) (|^) then the last formula (|^) 
obtains the very simple meaning: it uses the Fourier transform on [—tt, tt] 
to extract n-th term out of the generating function (I) for the polynomial 
sequence Pn{x). 


Remark 4.1 As we see that our formula (|^) can offer in conventional cal¬ 
culations the same (or slightly longer) procedure as the well-known umbral 
approach (|^). The situation could change if hypothetical quantum comput- 
will be able to calculate quantum propagators directly in parallel 


14 


ers 

processes. In this case the path integral formula should have advantages 
over the generating function (|]) where calculations must be done in a row. 

Finally we consider a realisation of the Schrodinger equation (|^. Ho¬ 
mogeneous propagator S{q,t) (U) are wave functions ([lOD itself. Thus they 
should satisfy to the Schrodinger like equation ([T^). For polynomial se¬ 
quences qn{x) this equation takes the form 


d ^ 

= Hqjx), 


(26) 


where the operator H is of the pseudodifferential type P^ , PB 

Han = TP^^ h{p)J^^^P an 


^-mpi 


— TV 
CO 


k=0 


h{p)J2<^ne^''^dp. 

^ ^ dk^n—ki 


fc=0 


(27) 
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where and the Fourier transform in the indicated variables. The 

function h{p) above is the Hamiltonian from the path integral and 

TT 

= e-^P-h{p)dp. 

— TT 

It turns to be just a convolution on N because Hamiltonian h{p) depends only 
fromp and is independent from k. The equation P7|) express the property 
of tokens to intertwine shift invariant operators between two cancellative 
semigroups. In the present case the operators are the derivative with respect 
X (i.e. the convolution with S'{x) on M) and the convolution with on M. 

Let M be the operator on sequences M : {a„,} h->■ {a„,/n!}. Then qn{x) = 
Mpn{x) and because M and ^ commute we obtain for a polynomial sequence 
Pn{x) of binomial type the equation 

^Pn{x) = M-^HMpn{x) ^pn(3:) =n\'^hk • (28) 

Of course the above equation follows from the differentiation of identity (|l|) 
with respect to y at the point p = 0. While this formula and the next 
example are rather elementary we emphasise their relations to the quantum 
mechanical framework. 




Example 4.2 The polynomial sequences from Examples |1.1| and |3.1| satisfy 
to the following realizations of equation 


1. The power monomials Pn{x) = x" satisfy to 

Pn-k{x) 


P'n{.x) = n\^5ki 


n\ 


k=0 


{n — k)\ {n — 1)! 


Pn-l{x) = npn-l{x). 


2 . 


3. 


The rising factorial sequence Pn{x) = 
to 

OO 

Pn{.x) =n\Y^ 
k=0 

The falling factorial sequence Pn{x) - 
to 


: x{x + 1) ■ ■ ■ (x + n — 1) satisfies 

1 Pn-k{x) 
k {n — k)\ 

= x(x — 1) ■ • • (x — n + 1) satisfies 


P'ni.x) 


I ( —1)^ Pn-k{x) 

k {n-k)V 
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4. The Abel polynomials An(x) = x{x — an)^ ^ satisfy to 

Pn-k{x) 


P: 


nix) = n! ^ 


k\ (n-k)\ 

k=0 ^ ' fc=0 


£ (t) ^Pn-kix) 

ju_n V / 


5. The Laguerre polynomials satisfy to 

CO /X 

/ / X I Pn-kix) 
fc =0 ^ ' 


Remark 4.3 In onr consideration we nsed primary the dehnition of to¬ 


ken (JT^) and the Pontrjagin dnality in ([T8|) -([I^ ) . Thns formnla has 
the same form for many other tokens [^, which we shall not consider here 
however. It follows from general properties of path integrals that for any 
admissible fnnction h{p) formnla gives a token qnix) and (P3|) prodnce 
a polynomial seqnence Pnix) of binomial type. 

We obtain fnrther generalisation if drop the homogeneity assnmption (|l^. 
Let combinatorial qnantities pn,kix,y) where n, k E N and x, y E M. he 
polynomials in the variable x of the degree n and in the variable y of the 
degree k. We assnme that they satisfy the identities 


Pn,kix,y) = ^Pn,mix,y)Pm,kiy,z), 


m=0 

for any fixed y snch that x < y < z. Obvionsly this is an inhomogeneons 
version of the identities (|^) and ([T^). The above qnantnm mechanical deriva¬ 
tion of path integral (^) can be carried (nnder some necessary assnmptions) 
also in this case. The main difference in the resnlting formnlae is that the 
Hamiltonian h{p, k) will depend from both variables p E [—vr, vr] and k E N. 
Conseqnently the Scrodinger eqnation (P^D will not be any more a trivial 
convolntion over N. We left consideration of this case for some fnrther pa¬ 
pers. 

In the conclusion: the main resnlt of this paper is yet another illnstration 
to the old observation that mathematics is indivisible. 
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